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I.     INTRODUCTION 

There  are  two  main  regions  of  deep  water  formation,   near-continent 
(Antarctic)  type  and  open  ocean  formation,   as  discussed  in  Killworth  (1979). 
The  Greenland  Sea  is  the  site  of  formation  of  the  largest  volume  of  deep  water 
in  the  North  Atlantic  (Helland-Hansen  and  Nansen,  1909).    Figures  (1)  and  (2) 
show  the  Greenland  Sea  area  of  interest  (Koltermann  and  Luthje,1989)  with  the 
significant  topographic  features  and  the  schematic  circulation. 

This  paper  will  examine  the  open  ocean  deep  water  formation  process  of 
the  Greenland  Sea.  Of  particular  interest  is  how  open-ocean  deep  water 
formation  is  influenced  by  the  depth-dependence  of  the  thermal  expansion 
coefficient  (<x).   It  is  hypothesized  that  free  convection  associated  with  parcel 
instability  is  possible  because  of  the  increase  in  the  thermal  expansion 
coefficient  with  pressure  in  conjunction  with  requisite  ambient  temperature  and 
salinity  profiles.  The  thermal  expansion  coefficient  increases  with  both  pressure 
and  temperature.  The  pressure  dependence  is  more  significant  for  cold  polar 
regions  than  it  is  for  the  deep  convection  of  the  Mediterranean,  as  shown  by 
Garwood  (1991).  The  dependence  of  a  on  pressure,  or  equivalents  depth,  is  a 
nonlinearity  in  the  equation  of  state  and  contributes  to  a  reduction  in  hydrostatic 
stability  (Carmack  and  Aagaard,  1973;  McDougall,  1987).   This  reduction  in 
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Figure  1 .   Greenland  Sea  Area  Chart  Area  of  interest  for  conditional 
instabilities. 
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Figure  2.  Greenland  Sea  Topography  and  Circulation:  topographic  features 
and  schematic  of  circulation. 


stability  only  occurs  when  the  gradient  of  temperature  is  negative  for  z  positive 
upwards. 

Parcel  instability  may  occur  if  a  water  parcel  that  is  cold  and  fresh  is 
displaced  downward  into  water  that  is  warmer  and  more  saline.   This  greater 
depth  causes  increased  a  for  a  surface  parcel  that  is  displaced  downward.   For 
a  sufficiently  large  change  in  a,  the  parcel  density  may  be  greater  than  the 
density  of  the  surrounding  water.   This  parcel  would  then  be  unstable  and 
would  be  accelerated  downward  by  gravity,  i.e.  a  parcel  instability  would  occur. 

The  traditional  explanation  of  deep-water  formation  (Nansen,  1906; 
Mosby,  1959)  focused  on  winter  cooling  of  surface  water  with  the  static  stability 
reduced  by  the  cyclonic  circulation  of  the  Greenland  Sea  gyre.   Killworth  (1979) 
hypothesized  that  with  progression  of  winter  cooling  the  depth  of  the 
overturning  water  column  increases  until  top-to-bottom  homogeneity  is 
achieved.   However,  there  have  been  no  observations  to  date  that  clearly 
support  vertical  homogeneity  over  the  whole  water  column. 

Carmack  and  Aagaard  (1973)  proposed  double  diffusion  as  the  cause  of 
the  deep  water  formation.   However,  Clarke  et  al.  (1990)  showed  that  the 
double  diffusion  process  can  not  account  for  the  observed  increases  in  bottom 
water  oxygen  concentration  over  expected  values  for  vertical  mixing  driven  by 
double  diffusion. 


Cabbeling  (Foster,  1972)  and  lateral  double  diffusion  have  been  proposed 
as  mechanisms  leading  to  deep  water  formation,  but  neither  of  these  processes 
are  evident  in  the  chimney  data  presented  by  Scott  and  Killworth  (1991) 
(Figure  3).   McDougall  (1984)  showed  how  diffusion  along  isopycnal  surfaces 
may  lead  to  water  mass  formation  by  both  cabbeling  and  thermobaricity.   This 
study  will  not  focus  on  isopycnal  diffusion,  but  on  surface-driven  buoyancy  flux. 

The  Greenland  Sea  stations  of  interest  are  located  in  the  center  of  the 
Greenland  Gyre  (GG)  north  of  the  Jan  Mayen  Current  (JMC).  Although 
Killworth  (1979)  believed  that  the  preconditioning  phase  does  not  occur  in  this 
region,  more  recently  Clarke  et  al.  (1990)  pointed  out  that  parts  of  the  Nordic 
Seas  are  considered  more  susceptible  to  deep  convection  than  are  other  areas. 
The  Greenland  Sea  region  circulates  cyclonically  and  is  weakly  stratified. 
These  conditions  permit  baroclinic  instability  (Killworth,  1979),  and  surface 
cooling  will  induce  vertical  mixing  by  creating  cooler  dense  water  over  less 
dense  water.   Killworth  (1979)  hypothesized  that  the  violent  mixing  initially 
occurs  during  winter  and  only  in  the  top  layers.   When  overturning  through  the 
water  column  results,  the  entire  water  column  has  characteristics  that  closely 
resemble  those  of  the  deep  water.  Aagaard  and  Carmack  (1989)  found  that  the 
deep  waters  in  the  Greenland  Seas  are  about  the  freshest  of  all  waters  found  in 
the  arctic  and  subartic.   They  proposed  that  the  surface  freshwater  plays  a 
larger  role  than  the  intermediate  water  in  the  formation  process.   Gascard 
(1990)  stated  that  deep  convection  which  leads  to  deep  water  formation  is  a 
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Figure  3.  Thermistor  Chain  Data:  Thermistor  chain  data  from  the  Greenland 
Sea  representing  a  chimney  event. 


surface-controlled  process  involving  intermediate  and  surface  waters.  The 
short  duration  of  about  one  week  for  chimney  events  suggested  by  Rudeis 
(1989)  may  explain  the  paucity  of  observations  of  convection  that  reach  the 
bottom  in  the  Greenland  Sea.   Lazier  (1973)  showed  an  isolated  example  of 
deep  convection  in  the  Labrador  Sea  reaching  approximately  1500  m,  but   it 
disappeared  within  2-3  days. 

Including  the  increase  in  the  thermal  expansion  coefficient  with  pressure, 
an  entrainment  model  has  been  developed  by  Garwood  (1991)  that  includes 
this  nonlinearity  in  the  equation  of  state.  This  model  predicts  greater  vertical 
mixing  and  penetrative  convection  than  do  earlier  models  because  of  the 
nonlinear  buoyancy  flux  enhancement.   In  the  study  here,  this  process  will  be 
examined  with  a  conventional  one-dimensional  model  applied  to  deep  oceanic 
mixing.  The  investigation  will  show  an  increase  in  vertical  turbulent  kinetic 
energy  that  may  enable  deep  penetrative  convection  and  the  hypothesized 
parcel  instabilities  leading  to  formation  of  deep  water. 

Before  applying  the  new  mixed  layer  theory,  however,  a  first  step  in 
Chapter  II  was  to  develop  a  "neutral  parcel"  model  that   predicts  a  neutrally- 
stable  temperature  profile  assuming  a  typical  two-layer  salinity  profile.  This 
profile  is  neutral  to  parcel  instabilities,  i.e.  an  isolated  parcel  of  water  from  the 
mixed  layer  would  have  a  neutral  buoyancy  when  displaced  vertically  anywhere 
in  the  water  column.  The  surface  temperature,  surface  salinity,  mixed  layer 
depth,  and  the  change  in  salinity  at  the  bottom  of  the  mixed  layer  are  specified. 


Then  the  neutral  parcel  model  diagnoses  the  change  in  the  temperature  at  the 
mixed  layer  bottom  and  the  change  of  temperature  with  depth  below  the  mixed 
layer.   Thus  for  a  given  two-layer  salinity  profile  a  temperature  profile  that  is 
neutrally  stable  for  displaced  mixed  layer  parcels  is  prescribed. 

The  next  step  after  development  of  the  neutral  parcel  model  is  to  study 
the  effects  of  the  depth  dependent  thermal  expansion  coefficient  on  mixed  layer 
dynamics.  A  calibrated  Kraus  and  Turner  (1967)  model  with  the  depth 
dependent  thermal  expansion  coefficient  is  applied  to  forced  and  free 
convection  to  demonstrate  the  effect  on  mixed  layer  dynamics. 

Chapter  III  will  expand  in  scope  from  the  mixed  layer  to  the  entire  water 
column.   Parcel  instabilities  in  the  water  column  are  examined  by  moving  a 
parcel  from  the  surface  or  from  within  the  water  column  over  the  entire  depth  of 
the  water  column.  Again,  the  thermal  expansion  coefficient  is  considered  to  be 
a  function  of  depth  and  temperatures  will  be  adjusted  for  heating  due  to 
compression  by  the  water  column.   Surface  parcels  moved  in  the  water  column 
will  then  be  tested  for  sensitivity  to  changes  in  initial  temperature  down  to 
freezing  and  then  increases  in  the  initial  salinity  due  to  brine  injection  from  ice 
formation. 

Chapter  IV  explores  regional  climatology  to  potential  parcel  instability, 
determining  regions  of  lower  stability.  Regional  stability  is  tested  by  examining 
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the  change  in  buoyancy  between  a  surface  parcel  and  the  in  situ  environment 
at  500  m  using  climatologicai  temperature  and  salinity  profiles. 


II.     NEUTRAL  PARCEL  STABIUTY 

A.      DESCRIPTION 

The  thermal  expansion  coefficient,  a,  is  defined  as 


p  ae 


(i) 


a«o0+ap  (2) 

where  p  is  the  relative  density  and  0  is  the  potential  temperature.   The  thermal 
expansion  coefficient  increases  with   pressure  (p)  as  a  nearly  linear  function  of 
pressure  (Garwood,  1991).   The   accurate   calculation  of  a  throughout  the 
water  column  requires  knowledge  of  the  ambient  salinity  and  the  ambient 
potential  temperature  versus  pressure  throughout  the  water  column.  The 
surface  value  of  the  thermal  expansion  coefficient  is  determined  by 


i  dp,      -p(s0,e0+Ae,p)+P(,s0,e0-Ae,p) 


(3) 


The  values  for  the  density  are  found  using  the  International  Equation  of  State 
for  seawater  (Millero  and  Poisson,  1981).   The  neutral-parcel  model  uses  a 
fixed  value  of  a^0.01  C.   The  centered  finite  differencing  scheme  is  used  along 
with  double  precision  due  to  the  small  differences  in  density  value.  An 
equation  similar  to  (3)  is  used  to  calculate  alpha  at  the  bottom  of  the  water 
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column  with  an  adjustment  to  the  potential  temperature  for  the  adiabatic 
temperature  increase  at  the  bottom  of  the  water  column. 

The  adjustment  for  the  adiabatic  temperature  rise  is  computed  using 
Bry den's  (1973)  algorithm  for  adiabatic  lapse  rate  as  a  function  of  salinity, 
pressure,  and  temperature.   By  prescribing  the  surface  potential  temperature, 
the  salinity  and  the  pressure  at  the  base  of  the  water  column,  the  increase  in 
potential  temperature  is  determined  using  Newton's  convergence  technique. 
The  temperature  increase  due  to  compression  is  added  to  the  surface  potential 
temperature  to  give  the  corrected  potential  temperature  at  the  base  of  the  water 
column.  With  the  calculated  temperature  and  given  salinity  and  pressure  at  the 
base  of  the  water  column,  the  bottom  value  of  alpha  (abottom)  is  calculated  using 
a  form  of  (3).   Equation  (3)  is  modified  for  otbottom  by  replacing  the  surface 
salinity  with  the  bottom  values  of  salinity,  pressure,  density  and  corrected 
potential  temperature.   Using  alpha  from  the  top  and  the  bottom  of  the  water 
column,  the  value  for  (a)  in  (2)  is  determined  using  (4).  Knowing  the  constants 


a_<bottoni)-a0  (4) 

botom 


of  (2),  the  change  in  the  value  of  alpha  from  the  top  to  the  bottom  of  the 
column  is  specified  as  a  linear  function  of  depth. 

The  expansion  coefficient  for  salinity  is  required  for  the  neutral-parcel 
model  to  evaluate  the  change  in  salinity  at  the  mixed  layer  interface  and  in 
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evaluating  the  temperature  values  in  the  modeled  "neutral  parcel"  water  column. 
The  value  of  p  is  determined  using  (5).   The  expansion  coefficient  for  salinity,  (3, 

p=JL  3^0799  (5) 

p  d$        p 

is  assumed  constant  throughout  the  water  column.  The  value  of  0.799  psu  1  is 
the  average  of  p  values  over  a  400  bar  water  column  (Gill,  1982). 

With  the  above  information,  the  neutral-parcel  model  can  now  be  used  to 
prescribe  the  values  for  the  potential  temperature  profile  as  a  function  of  depth 
or  pressure.   The  physical  meaning  of  the  neutral-parcel  model  is  that  the 
profile  is  stable  everywhere  for  small  vertical  displacements,  with  Brunt-Vaisala 
frequency  (N2>0).  However,  a  surface  parcel  that  is  moved  downward  will  be 
perfectly  neutral,  remaining  wherever  it  is  placed  with  no  restoring  force. 
For  a  neutral  buoyancy  profile,   (6)  is  used  to  solve  for  the  potential 
temperature  as  a  function  of  depth,  (7). 

P(2)=P[1-«W[e-0j-p[5-5j]  (6) 

e«=e0+-^  (7) 

Figure  (4)  shows  profile  data  for  station  82  of  the  Marginal  Ice  Zone 
Experiment  (MIZEX-87)  (Quadfasel  et  al.,  1988).   The  data  from  MIZEX-87  was 
collected  during  the  RV  "Valdivia"  cruise  54  in  March-April  1988  in  the 
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Figure  4.   Neutral-Parcel  Model  Run  on  Station  82:  The  top  two  panels  are  the 
T-S  profiles  of  Station  82.   The  bottom  panels  show  the  neutral-parcel  model 


run. 
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Greenland  Sea  and  the  Fram  strait.  One  objective  of  this  cruise  was  to  find 
and  map  an  active  convective  element.   Station  82  has  representative 
temperature  and  salinity  profiles  for  the  center  of  the  Greenland  Gyre.  In 
Figure  (4),  the  lower  right  plot  (DIFFERENCE)  represents  the  difference 
between  station  82's  temperature  profile  and  the  neutral  parcel  temperature 
profile.  This  "difference"  demonstrates  how  close  station  82's  temperature 
profile  is  to  a  neutral  parcel  temperature  profile.  This  also  demonstrates  that 
the  center  of  the  Greenland  Gyre  is  close  to  neutral  stability. 

B.      EFFECT  OF  ot(p)  ON  MIXED  LAYER  DYNAMICS 

To  demonstrate,  the  effect  of  adding  a(p)  to  a  calibrated  Kraus  and 
Turner  (1967)  model,  three  cases  are  examined:  forced  convection,  free 
convection  and  a  combination  of  forced  and  free  convection.   Using  the  Kraus 
and  Turner  model  (8)  and  (9)  and  the  constants  listed  in  Table  I  determined 


{-figAS+*gh*  7)^=0,^-^^-%) 
dt  pCp 

dS_{F+E-Pis ar_ 

at'       h  h 


(8) 


ATdh 
BT_   Qp  _        dt  O) 

dt     pCf/l        h 
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TABLE   1.      CONSTANTS    IN  THE   PARCEL   INSTABILITY  MODEL 

I 

1 . 5    cm/s 

g 

983    cm/s2 

Ah0 

1800    cm 

AT0 

+    1.52    C 

ffo 

0.475x10-"    s'1 

a, 

0.265xl0"7  cm^s"1 

Qo 

-0.005   cal/cm 

P 

1   g/cm3 

CP 

1    cal/g/K 
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from  the  NPS  mixed  layer  model  (Garwood,  1977)  for  arctic  conditions,  the 
results  of  forced  and  free  convection  for  both  a  constant  value  of  a,  and   a  as 
a  function   of  depth   will  be  solved  numerically  and  compared.   In  (8),  F  is  the 
rate  of  water  freezing,  E  is  the  rate  of  evaporation  and  P  is  the  rate  of 
precipitation.   In  (8)  and  (9),  h  defines  the  mixed  layer  depth,  AT  is  the  change 
in  temperature  at  the  mixed  layer,  g  is  gravity,  and  Q0/pcp  is  the  downward 
surface  temperature  flux.   Constants  c,  and  c2  are  dimensionless  model 
constants. 

The  cooling  case  (Q0<  0)  with  continued  deepening  of  the  mixed  layer 
over  time  will  be  considered  first.  The  left  hand  side  of  (8)  is  the  rate  turbulent 
kinetic  energy  (TKE)  is  converted  to  potential  energy  (PE)  by  entrainment  (layer 
deepening).   The  first  term  on  the  right  hand  side  of  (8)  is  the  generation  rate  of 
the  TKE  by  wind  stirring  minus  viscous  dissipation.   The  second  term  is 
damping  (Q0>  0)  of  TKE  by  surface  heating  or  free  convection  production 
(Q0<  0)  of  TKE  by  surface  cooling.  The  left  hand  side  of  (9)  is  the  rate  of 
change  of  the  temperature  at  the  bottom  of  mixed  layer.   The  first  term  on  the 
right  hand  side  of  (9)  is  the  surface  heating  and  the  second  term  is  the 
entrainment  heat  flux,  with  a  mixed  layer  temperature  decrease(increase)  if  AT 
is  positive(negative). 
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For  case  1 ,  forced  convection  (Q0=  0),  the  wind  forcing  is  the  only  source 
of  turbulent  kinetic  energy.   For  this  case  (8)  and  (9)  simplify  to  (10)  and  (11). 


(10) 


dt 

A  73/7 
h   dt 

W? 

dh_         C|"»  (11) 

dt    agh±T-$gAS 


Note  that  hAT  equals  a  constant: 


A7=7-7{Z=-/7) 
3A7_37_37(/Tv 

dt    dt   dr   ' 

4^(-/7)=0 

dr   ' 

37_    A73/J 

dt"  h  dt 

h*±I+ A7^?=0 
dt  dt 

dr     ' 


For  a  constant  over  depth  and  a  positive(negative)  AT,  integrating  (11)  from  hc 
to  h  will  result  in  a  linearly  increasing(decreasing)  mixed  layer  depth  as  a 
function  of  time. 
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With  the  same  initial  conditions,  (11)  becomes  (12)  if  a  is  depth 
dependent.   Integrating  from  h0  to  h  results  in  (13).   Figure  (5)  shows  the 
deepening(shallowing)  of  the  mixed  layer  depth  for  constant  a  and  depth 
dependent  a  versions  of  case  1 .  The  mixed  layer  deepening(shallowing)  is 
slowed  over  a  one  day  period  when  a  is  depth  dependent.   Furthermore,  the 

dh_    C,(/  1  (12) 

dt   gh^T0\0-a,z/ 

h-h=-^- 1—  (13) 

0  9hJLT0{*satH 

depth  dependent  a  causes  a  nonlinear  change  in  depth  of  the  mixed  layer  as 
noted  by  the  slight  curvature  of  the  dotted  line  for  positive  AT  and  the  dashed 
line  for  negative  AT  in  Figure  (4). 

Case  2,  free  convection  (u.=  0),  the  condition  of  surface  cooling  (Q0<  0) 
for  a  constant  a  and  a(z)  are  compared.   Using  the  Kraus  and  Turner  model  on 
the  free  convection  case  results  in  (14).   Notice  in  (14)  that  a  appears  in  both 


3/7  = Pfy_ 

dt        aghAT 
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Figure  5.   Case  1:  Forced  convection  mixed  layer  deepening(shallowing), 
constant  a  is  the  dot-dashed(solid)  line  and  a(z)  is  the  dotted(dashed)  line. 
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the  numerator  and  the  denominator.   It  is  not  eliminated  because  a  now  has 
different  rates  at  the  top  and  bottom  of  the  mixed  layer.  The  first  term  on  the 
right  hand  side  of  (15)  is  the  turbulent  heat  flux  at  the  surface,  TW  [.   The 
second  term  on  the  right  hand  side  of  (15)  is  turbulent  heat  flux  at  the  bottom  of 
the  mixed  layer,  Tw'  [,. 

/^A-AT^  (15) 

dt   pcp       dt 

When  ocwas  assumed  constant,  the   a  values  in  the  numerator  and 
denominator  can  be  divided  out  which  reduces  (14)  to 


dt      ^  PC. 


And  recall  from  earlier  that  at  z=-h,  3  A  T_  ST 

~3T~~3f 
Then, 

dt  dt      pcp 

dt  pCD 


Now  integrating  from  h0  to  h  results  in  -CJQ 

hA  T-h,A  T0= 1 

Solving  (16)  numerically  gives  solutions  in  Figure  (6). 
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For  the  depth  dependent  a  case,  substitute  asa0-a,z  into  (14).   The  value 
of  a  in  the  numerator  will  be  a0  since  it  is  the  surface  flux  component.   The 
denominator  a  value  becomes  ot0+  a,  h(t)   at  the  bottom  of  the  mixed  layer 
(z  =  -h).   Using  these  substitutions  and  solving  at  z  =  -h  yields, 


|(/,A7)^>(— £— ) 

dt  pcp    ao+a,/<0 

Integrating  from  h0  to  h  gives  (17).   Figure  (6)  shows  that  the  mixed  layer  is 


hAT-h^T.^&f—ll—dt  (17) 


deepening(shallowing);  however,  the  rate  of  deepening(shallowing)  is  much 
slower  than  in  case  1 .   Note  that  the  deepening(shallowing)  is  much  slower 
than  that  in  case  1 .  For  the  a(z)  case,  the  magnitude  of  reduction  in  the  rate  of 
mixed  layer  deepening(shallowing)  is  smaller  than  that  in  the  forced  convection 
case.  Thus  the  results  show  that  as  in  the  forced  convection  case,  the  slowing 
of  the  mixed  layer  deepening(shallowing)  is  caused  by  the  increase  in  a  with 
time. 
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Figure  6.   Case  2:   Free  convection  mixed  layer  deepenlng(shallowlng), 
constant  a  is  the  dot-dashed  (dotted)  line  and  a(z)  is  the  dashed(solid)  line. 
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Cases  1  and  2  were  a  simple  analysis  of  how  the  Kraus  an  Turner  (1967) 
model  was  effected  by  simply  plugging  in  a(z).   Case  3  will  examine  how  the 
Kraus  an  Turner  (1967)  model  is  effected  when  a(z)  is  considered  in  the 
derivation  of  (8)  and  (9).   For  case  3,  the  combined  effects  of  both  forced  and 
free  convection     (Q0(  0  and  u.  *  0)  are  examined.   The  first  term  on  the  right 
hand  side  of  (8),  generation  rate  of  wind  stirring  minus  viscous  dissipation,  can 
be  rewritten  in  terms  of  turbulent  heat  flux  as  shown  in  (18).  The  top  line  of  (8) 
can  be  rewritten  as  (19)  using  form  shown  in  (18),  the  generation  rate  of  TKE 
balance  by  the  vertically  integrated  buoyant  dampening. 

s  1  r°    ^7/w-  (18) 


<VJ?=lf°h*gT'w>dz 


*A1h£4/>rvto  (19) 


The  turbulent  heat  flux  is  a  function  of  the  depth  and  is  written  in  the  form 
Multiply  the  turbulent  heat  flux  by  ag,  substitute  in  the  expressions  for  the 


T'w'{% = T  V(0)(1  +4)  -  T  V(-ty(4)  (20) 

h  h 

surface  turbulent  heat  flux  and  the  turbulent  heat  flux  at  the  bottom  of  the 
mixed  layer,  and  integrate  over  the  depth  of  the  mixed  layer.  These  operations 
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result  in  (21).  The  left  hand  side  of  (21)  is  the  integral  average  of  the  buoyancy 
flux  over  the  water  column.  The  first  term  on  the  right  hand  side  is  the 
buoyancy  flux  component  due  to  the  net  downward  temperature  flux.  The 
second  term  is  the  buoyancy  flux  component  due  to  the  bottom  heat  flux. 


/>^'*-/>1+|M-#)*+/>fi(|)(A7f*       «21» 


For  constant  a,  we  can  reduce  (21)  by  dividing  out  the  gravity  and  solving 
the  integral  to  give 


If  °arV(fe-[-^-A^  (22) 

hi-h  l    ocn         bV  2 


PC, 

For  the  depth  dependent  a,  we  will  substitute  into  (21)  the  form   a  =  a0  -  a,z 
from  (2): 

f  °agTfw'dz 

j  -h 

-/>,--.4(i+f>(-^)* 

♦/ >«o-«.4(f)(A  T^)dz  | 

Expanding  the  integral  gives  (23). 
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Simplify  (23)  results  in, 


^>*W*^Sf>-MF»*       ,23) 


Evaluating  this  expression  gives  (24);  the  buoyancy  flux  components  due  to  net 
downward  temperature  flux  and  bottom  heat  flux. 

■*£*M*-*IfM*  (24) 

In  comparing  (22)  and  (24),  notice  that  the  depth  dependent  a  causes  an 
increased  weighting  to  both  the  temperature  and  bottom  heat  flux  components. 
Another  difference  is  the  weighting  between  the  two  fluxes.   The  bottom  heat 
flux  has  a  larger  weight  which  makes  it  a  larger  reduction(addition)  in  the 
buoyancy  flux  due  to  the  net  downward  temperature  flux,  depending  on  whether 
AT  is  positive  or  negative.   This  shows  that  the  depth  dependent  a  causes  the 
mixed  layer  to  deepen(shallow)  at  a  slightly  slower  rate  for  a  positive(negative) 
AT  than  what  is  found  for  constant  a. 
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III.     TEST  OF  PARCEL  INSTABILITY  THEORY 

A.      WATER  COLUMN  INSTABILITIES 

The  response  of  the  mixed  layer  to  the  nonlinear  effects  of  ot(z)  was 
investigated  in  the  previous  section,  but  the  possibility  of  instability  below  the 
mixed  layer  needs  to  be  considered.   Figure  (3)  depicts  Scott  and  Killworth's 
(1991)  thermistor  chain  data  from  a  track  between  Greenland  and  Iceland.   The 
temperature  contours  show  the  existence  of  two  chimneys  below  the  surface. 
The  chimneys  represent  events  that  are  clearly  below  the  surface  boundary 
layer.  Are  there  features  analogous  to  cumulus  cloud  formation  in  the 
atmosphere?   Riehl  (1979)  describes  how  atmospheric  deep  convection  in  the 
tropics  requires  only  a  small  percentage  of  the  region  to  be  conditionally 
unstable.  Is  this  the  same  situation  seen  in  deep  chimney  oceanic  convection? 
Scott  and  Killworth  (1991)  hypothesized  that  the  chimneys  in  Figure  (3)  are 
decoupled  from  the  surface.  A  chimney  could  be  coupled  with  the  surface  if  it 
were  created  by  a  small  instability  in  a  region  of  conditional  instability.  A  small 
instability  from  the  surface  would  not  be  very  evident  once  the  chimney  event 
started.   Thus  what  looks  decoupled  after  the  chimney  is  developed  may  not 
show  how  the  surface  initiated  the  event. 

To  determine  how  an  instability  in  a  region  of  conditional  instability  relates 
to  our  neutral-parcel  instabilities,  the  model  is  applied  to  station  data  in  the 
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Greenland  gyre.   By  displacing  a  surface  parcel  downward  in  the  water  column, 
the  possibility  of  parcel  instability  is  explored.   The  surface  parcel  may  be 
cooled  and  its  salinity  may  also  be  changed  by  freezing.   The  effect  on  stability 
for  such  alleged  surface  parcels  is  also  examined.   This  procedure 
demonstrates  what  changes  in  a  surface  parcel  can  cause  instability  or  make 
the  parcel  neutrally  stable/unstable  and  therefore  favorable  for  a  convective 
event. 

B.      MODEL  DESCRIPTION 

To  investigate  the  possibility  of  parcel  instabilities  using  in  situ  data, 
temperature  and  salinity  values  from  station  44  of  MIZEX  (Quadfasel  et  al., 
1988)  are  applied  to  the  parcel  instability  model.   The  temperature  and  salinity 
data  is  used  in  the  Millero  and  Poisson  (1981)  equation  of  state  to  obtain 
density  values  throughout  the  water  column.   The  temperature  is  adjusted  for 
the  heat  of  compression  with  increasing  depth  by  using  Bryden's  (1 973) 
algorithm  to  obtain  potential  temperature  (degrees  C)  as  a  function  of  salinity 
(psu),  temperature  (C)  and  pressure  (bars). 

The  in  situ  density   is  calculated  at  each  level  of  the  temperature  and 
salinity  profile.  Then  moving  the  surface  parcel  down  the  water  column  profile, 
the  temperature  is  adjusted  for  compression  at  each  new  pressure  level  using 
the  Bryden  algorithm.  With  the  compression-adjusted  surface  temperature  and 
surface  salinity,  the  density  is  calculated  for  the  surface  parcel  at  each  pressure 
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level.   Each  level  has  an  in  situ  density  and  a  density  for  the  surface  parcel 
moved  to  that  level.   The  change  in  density  is  found  at  each  level  by 

When  Aps  is  negative,  the  surface  parcel  would  be  unstable  (more  dense) 
compared  to  its  surroundings  and  would  accelerate  downward  by  gravity  until  it 
reached  a  depth  with  equal  density  or  until  it  encountered  the  bottom. 

Station  44  will  be  examined  to  determine  if  the  situation  described  above 
ever  occurs.   The  results  are  seen  in  Figure  (7),  a  plot  of  depth  versus  Ap. 
The  solid  line  shows  the  non-compression  adjusted  temperature  values  for  the 
change  in  density  at  each  level  as  depth  increases.  The  dashed  line  shows  the 
adiabatic  temperature  adjusted  surface  parcel  changes  in  density  as  depth 
increases.  The  profiles  do  not  exhibit  any  parcel  instabilities  for  the  displaced 
surface  parcel  since  the  values  of  Ap  never  become  negative.  Figure  (7)  also 
shows  that  the  temperature  adjustment  is  a  stabilizing  effect  for  the  profile;  the 
adiabatic  temperature  corrected  Ap  curve  is  more  positive  through  the  water 
column. 

1 .       SENSITIVITY  TO  DEPTH 

The  procedure  described  is  then  modified  to  begin  with  a  parcel  at 
the  bottom  of  the  mixed  layer  vice  a  surface  parcel.  No  parcel  instabilities  were 
noted.  Other  levels  are  also  examined,  and  no  parcel  instabilities  were 
discovered.   Figure  (8)  is  an  example  of  this  procedure.   The  dot-dash  line 
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Figure  7.   Change  In  density  curves  for  Station  44:  The   dotted  curve 
represents  adiabatic  temperature  adjustment.   The  solid  curve  is  not  adiabatic 
temperature  adjusted. 
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Figure  8.   Change  in  density  curve:  Testing  for  depth  sensitivity,  solid  being 
non-adiabatic,  dotted  being  adiabaticly  adjusted  and  dot-dash  is  a  500  m  parcel 
moved  in  the  water  column. 
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shows  a  parcel  from  500  meters  depth  moved  to  the  surface  and  then  moved 
down  the  water  column.   Below  500  meters  the  500  meter  parcel  never  became 
unstable. 

It  is  worth  noting  that  Gascard  (1991)  stated  that  chimney  formation 
only  occurs  were  the  overall  stratification  can  be  destabilized  and  destroyed 
from  the  interior.   He  further  stated  that  vertical  mixing  by  turbulent  entrainment 
triggered  from  the  surface  is  ineffective.  Figure  (8)  shows  an  example  of  a  test 
for  destabilization  of  the  stratification  from  the  interior,  i.e.  500  meters  in  this 
example.   However,  looking  at  the  above  parcel  instability  model  for  the  data  in 
the  Greenland  Sea  given  by  Gascard  (1991)  shows  no  parcel  instabilities  using 
water  from  intermediate  levels.   Figure  (9)  demonstrates  how  a  1000  meter 
parcel  moved  in  the  water  column  did  not  cause  instability  either. 

2.       SENSITIVITY  TO  SURFACE  TEMPERATURE 

The  next  test  was  to  assume  that  the  surface  was  cooled  in  steps 
down  to  the  freezing  temperature  by  either  heat  loss  to  the  atmosphere  or  by 
advection.    By  using  Millero's  (1978)  equation  for  the  freezing  point  of 
seawater,  the  maximum  amount  of  cooling  that  could  be  imposed  on  the 
surface  is  calculated.  At  station  44,  the  surface  temperature  was  already  - 
1.797  C  and  could  only  be  cooled  by  0.1  C  to  reach  the  freezing  point. 

The  first  step  is  to  apply  the  parcel  instability  model  with  station  44's 
original  surface  temperature.   Each  consecutive  test  is  applied  to  the  same 
station  44  water  column  profile  but  with  the  surface  parcel  temperature 
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Figure  9.   Change  in  density:  Test  of  Gascard's  Greenland  Sea  data  for 
internally  caused  instabilities. 
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Figure  10.   Contour  of  Ap(0,Z):  Sensitivity  to  reducing  surface  temperature  to 
freezing. 
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reduced  by  0.02  C  at  the  start  of  each  iteration  until  the  freezing  temperature 
was  reached.   Figure  (10)  shows  the  contours  of  Ap(0,z).  Viewing  Figure  (10) 
from  right  to  left  shows  the  result  of  decreasing  station  44's  surface  parcel 
temperature  to  the  freezing  point. 

Figure  (11)  shows  an  expanded  view  of  the  upper  750  m  for  the 
same  contour  intervals.  The  only  instabilities  demonstrated  are  in  the  upper  50 
m  of  the  water  column,  the  mixed  layer.   One  interesting  feature  to  note  is  a 
zone  of  maximum  stability  between  500  m  and  1500  m  that  is  strongest  at  the 
original  surface  parcel  temperature  and  becomes  less  dominant  as  surface 
parcel  temperature  approaches  freezing.  This  can  is  evident  from  the  0.08  Ap 
contour  in  Figure  (11). 

3.       SENSITIVITY  TO  SURFACE  SALINITY 

With  no  parcel  instabilities  apparent  in  the  water  column  when  the 
temperature  was  reduced  to  the  freezing  point,  the  effects  of  surface  salinity 
changes  are  examined.    An  increase  in  salinity  is  associated  with  freezing  after 
the  surface  water  is  cooled  to  the  freezing  point.   Then  the  change  in  salinity 
(8S)  times  the  thickness  of  the  ice  (h)  is  proportional  to  the  change  in  the  ice 
thickness  (6h)  times  the  difference  in  surface  salinity  (S)  and  ice 
salinity  (S,), 

The  salinity  of  the  ice  is  approximated  with  Maykut's  (1985)   ice  salinity  data. 
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Figure  11.   Contour  Plot  of  Ap(0,Z):  Expanded  view  of  the  upper  750  meters  of 
sensitivity  to  reducing  surface  temperature  to  freezing. 
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Figure  12.   Contours  of  Ap(S,z):  Sensitivity  to  increasing  salinity  by  .02  psu 
with  the  surface  temperature  at  freezing. 
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An  ice  salinity  value  of  12  psu  for  first  year  ice  of  initial  thickness  at  5 
cm  or  less  will  give  approximately  0.1  psu  increase  in  salinity  for  1  cm  of  ice 
growth  at  the  surface  of  the  water  column.   This  does  not  account  for  the 
salinity  increases  due  to  mixed  layer  turbulent  mixing  or  advection  of  colder 
more  saline  water  into  the  surface  water.  Anderson's  (1961)  data  for  ice  growth 
rate  gives  a  10  cm/day  ice  growth  rate  for  a  storm  with  air  temperatures  at  -40 
C.  Thus  an  increase  of  0.1  psu  for  the  surface  water  is  reasonable  to  use  in 
the  parcel  instability  model. 

Solving  the  parcel  instability  model  for  the  surface  parcel  at  freezing 
temperature  and  increasing  the  salinity  by  0.1  psu  in  0.02  psu  increments  gives 
the  contour  plot  of  figure  (12).   The  graph  is  salinity  vs.  depth,  plotting  the 
change  in  density  between  the  surface  parcel  moved  to  depth  as  salinity 
increases  by  0.02  psu  increments,  Ap(salinity .depth).  As  the  salinity  is 
increased  from  its  original  value  at  the  surface  of  34.74  to  34.80  psu  and 
greater,  the  parcel  instability  (negative  values  of  Ap)  are  predicted.  Also,  values 
greater  than  34.84  psu  show  parcel  instability  to  the  bottom  of  the  water 
column.  The  second  to  the  last  iteration  (second  to  last  contour  on  the  right) 
shows  the  neutral  stability  curve  for  a  parcel  moved  from  the  surface  to  depth 
without  any  change  in  the  density  between  the  surface  parcel  at  the  freezing 
point  and  the  in  situ  water  density. 

Figure  (13)  shows  an  expanded  view  of  the  upper  750  m  for  the 
same  contour  intervals.   This  shows  how  the  mixed  layer  instability  grows  very 
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Figure  13.   Contour  Plot  of  Ap(S,z):  Expanded  view  of  upper  750  meters, 
sensitivity  to  increasing  salinity  with  surface  temperature  at  freezing. 
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quickly  from  the  onset  of  the  salinity  increase.   Examining  the  values  of 
instability  near  the  salinity  value  of  34.76  psu  shows  at  what  value  parcels 
should  "break  out"  of  the  mixed  layer  in  deep  convection  events.   Thus  the 
model  demonstrates  that  a  parcel  instability  can  occur  for  the  station  44  profile 
provided  that  the  surface  is  cooled  to  freezing  and  as  little  as  1  cm  of  ice  forms. 
The  salinity  injection  and  surface  temperature  changes  at  the  surface  can 
create  a  parcel  instability  that  can  cause  convection  to  the  bottom  of  the  water 
column. 
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IV.     CLIMATOLOGY 

A.  CLIMATOLOGY  AND  PARCEL  INSTABILITY 

Is  there  a  region  in  the  Greenland  sea  that  can  be  considered  preferential 
for  parcel  instability  events  leading  to  chimneys?  With  the  assumption  that 
parcel  instabilities  depend  on  a  depth  dependent  thermal  expansion  coefficient 
the  Levitus(1982)  Climatological  Atlas  of  the  World  Oceans  is  used  to 
investigate  the  Greenland  Sea.   The  Levitus  (1982)  data  set  contains 
temperature  and  salinity  at  standard  oceanographic  observation  levels  on  a  one 
degree  latitude  and  longitude  grid  for  the  world  oceans.  The  data  are  averages 
of  objectively  analyzed  station  data,  mechanical  bathythermograph  and 
expandable  bathythermograph  between  the  surface  and  a  maximum  depth  of 
5500  meters.  Annual  and  seasonal  data  were  provide,  but  only  the  winter 
season  data  is  used  for  the  parcel  instability  analysis  here. 

B.  CLIMATOLOGICAL  MODEL  TEST 

Again  using  the  depth  dependent  thermal  expansion  coefficient  and  the 
principles  of  the  parcel  instability  model,  the  temperature  and  salinity  values 
from  55  N  to  85  N  in  latitude  and  5  E  to  20  W  in  longitude  will  be  used  to 
determine  the  stability  of  the  region.   The  surface  values  of  temperature  and 
salinity  from  the  data  set  do  not  represent  a  very  realistic  mixed  layer  due  to 
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the  averaging  process,  but  it  is  assumed  they  provide  an  adequate  first  guess 
field.  A  depth  of  500  m  will  be  the  comparison  depth  for  the  instability  test. 
This  depth  is  used  to  help  cover  the  largest  area  of  the  data  set  without  being 
bottom  limited.  The  500  m  depth  is  also  a  reasonable  chimney  depth  evident  in 
various  data  sets. 

The  model  will  use  (25)  to  determine  the  buoyancy  jump  (Ab)  as  a 
function  of  the  500  m  depth  to  create  a  grid  of  the  change  in  buoyancy  ( Ab) 

A^fl^TWrW-P^-iW)  (25) 

values  to  contour.  The  values  subscripted  with  500  indicate  the  value  at  500  m 
and  the  subscript  surf  indicates  the  surface  value.  Again,  in  this  test  (3  will  be 
kept  constant. 

Figure  (14)  represents  the  Ab  values  plotted  on  the  latitude  and  longitude 
grid.   The  area  encircled  by  the  0.1  cm/s"  contour  represents  an  area  close  to 
neutral  stability.   Comparison  of  this  area  to  the  charts  in  Figures  (1)  and  (2) 
show  that  this  is  highly  correlated  to  the  central  region  of  the  Greenland  Gyre. 
The  areas  on  Figure  (14)  without  a  contour  value  represent  land  or  the  water 
depth  was  too  shallow  for  comparison  with  the  500  m  depth. 
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Figure  14.    Contoured  Plot  of  Change  in  Buoyancy:  Contours  of  the  Ab 
between  the  surface  values  and  a  depth  of  500  meters. 
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This  instability  test  supports  the  theory  that  the  Greenland  Gyre  is  the 
most  likely  candidate  for  chimney  formation  in  the  Greenland  Sea  region.   It 
helps  provide  a  bounds  for  further  data  collection  and  observations  of  chimney 
events. 
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V.     SUMMARY 

This  study  has  (1)  investigated  the  dynamics  of  a  neutral-parcel  model 
with  a  depth  varying  thermal  expansion  coefficient  (a),  and  (2)  tested  for  parcel 
instability  in  the  Greenland  Sea,  using  both  actual  profiles  and  climate-logical 
data. 

In  the  neutral-parcel  model,  the  effects  on  mixed  layer  dynamics  for  forced 
and  free  convection  and  the  combination  of  both  were  analyzed  for  what  added 
dynamical  effect  the  depth  dependent  a(z)  had  on  the  mixed  layer  in  each  test. 
The  depth  dependent  a(z)  changes  significantly  mixed  layer  deepening.  For  all 
three  cases  of  convection,  depth  dependent  a(z)  caused  the  mixed  layer  to 
deepen(shallow)  at  a  slower  rate  given  a  positive(negative)  AT  than  what  was 
found  for  a  constant  a. 

The  parcel  instability  model  using  a(z)  and  adiabatic  temperature 
increases  due  to  parcel  movement  down  the  water  column  was  tested  for 
sensitivity  to  parcels  displaced  from  interior  depths,  sensitivity  to  reducing  the 
surface  temperature  to  freezing  and  sensitivity  to  increases  in  the  surface 
salinity.  Gascard's  (1991)  destabilization  from  the  interior  could  not  be 
supported  using  the  parcel  instability  model,  demonstrated  by  the  test  for  depth 
sensitivity  using  data  from  Gascard  (1991).   Reducing  the  surface  temperature 
to  freezing  showed  the  onset  of  freezing  reduced  the  stability  and  brought  it 
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closer  to  neutral  stability.    Cooling  surface  temperature  to  freezing  did  not  by 
itself  lead  to  a  state  of  instability  for  most  of  the  profiles  tested.     When  the 
surface  was  at  freezing  and  enough  ice  formation  occurred,  then  the  increase 
in  salinity  resulted  in  the  possible  onset  of  parcel  instability. 

Using  the  depth  dependent  thermal  expansion  coefficient  and  adiabatic 
temperature  increases,  an  assessment  of  the  possible  instability  in  the 
climatological  data  of  the  Greenland  Sea  was  made.  The  area  of  least  stability 
for  this  test  was  the  center  of  the  Greenland  Gyre. 
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